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Abstract 



We investigate the influence of a time dependent, homogeneous electric field 
on scattering properties of non-interacting electrons in an arbitrary static 
potential. We develop a method to calculate the (Keldysh) Green's function in 
two complementary approaches. Starting from a plane wave basis, a formally 
exact solution is given in terms of the inverse of a matrix containing infinitely 
many 'photoblocks' which can be evaluated approximately by truncation. In 
the exact eigenstate basis of the scattering potential, we obtain a version of 
the Floquet state theory in the Green's functions language. The formalism 
is checked for cases such as a simple model of a double barrier in a strong 
electric field. Furthermore, an exact relation between the inelastic scattering 
rate due to the microwave and the AC conductivity of the system is derived 
which in particular holds near or at a metal-insulator transition in disordered 
systems. 

PACS numbers: 71.10 +x , 72.30 +q 
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I. INTRODUCTION 



The influence of a time-dependent electric field on properties of a quantum meclianical 
system is tlie topic of active recent researcfi. An essential point in these investigations is 
that the time-dependent field is not from the beginning treated in perturbation theory (e.g. 
in linear response), but is rather considered as inherent part of the system itself. By this, 
one has to deal with conditions of a non-equilibrium situation under which the quantity of 
interest, e.g. a tunnel current or the response to an additional DC field, has to be determined. 
There is a considerable variety of physical systems studied so far. Many works are related 
to tunneling, such as driven quartic double wells [|I|, driven quantum wells ^ or tunneling 
through laser irradiated wells 0. In interacting systems, investigations concentrated on 
tunneling through quantum dots in time- varying fields 

Another class of investigations are related to transport experiments as, e.g., in the quan- 
tum Hall regime [^|-|^, where a microwave gives rise to peculiar changes of the DC conduc- 
tivities. Calculations have been done on photo-induced transport through quantum point 
contacts §], or on microwave induced changes of the DC conductivity of a disordered system 
T0| . Recently, Holthaus and Hone studied the Wannier-Stark ladder in an ac-driven 



one-dimensional tight- binding system, where in presence of disorder-induced localization 
the localization length is controlled by the time-dependent field. 

In spite of this great variety of investigations, the main theoretical focus so far has been 
on the zero- or (quasi) one- dimensional case. It is therefore desirable to develop a calculation 
scheme that is applicable in higher dimensions and for arbitrary static potentials, in partic- 
ular in view of applications as the quantum Hall effect {d = 2) or the Anderson transition 
{d = 3). In this paper, we set up a general formalism to calculate the Green's function of 
non-interacting electrons moving in an arbitrary static potential under the influence of a 
time-dependent electric field with frequency Q. We use the Dyson equation to obtain the 
Keldysh Green's function in two different calculation schemes. First, we regard the static 
potential as a perturbation and include the electric field in the unperturbed Hamiltonian. 
We give a formally exact solution which can be evaluated approximately by inverting a trun- 
cated matrix containing a finite number of 'photo-blocks'. The advantage of this method 
is its exactness in the electric field; it furthermore sums up the static potential to infinite 
order and is perturbative in the higher Fourier components of the Green's function which 
correspond to the 'center-of-mass' time coordinate. Numerical examples of this 'truncation 
method' are presented for the trivial case of a constant potential, which allows comparison 
to an exact result, and a simple model for the transmission through a double barrier struc- 
ture in a strong electric field. There, our method qualitatively reproduces previous results 
by other authors |T^J§]. In particular, the examples show that already a small number of 
'photo-blocks' is sufficient to obtain convergence. 

Second, starting from the exact eigenstates of the static potential, again an exact formal 
solution is derived in which the Green's function is represented as the inverse of an infinite 
tridiagonal matrix, which is the Green's function analogue of the Floquet state Hamiltonian. 
This approach in particular is useful in situations where the static scattering problem is 
already solved (say by numerical diagonalization) and one is interested in the effect of an 
additional, time-dependent electric field. We calculate the electron self energy of which the 
diagonal part yields an inelastic scattering time tmw due to the influence of the electric 
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field. In the limit of hfl <^ Ep (where Ep is the Fermi level of the system without electric 
field) we obtain an exact relation between tmw a-nd the zero temperature AC-conductivity 
of the system without electric field. We calculate tmw for microwave radiation in the integer 
quantum Hall effect. From the fact that tmw alone is insufficient to interpret experimental 
data, we conclude that one rather has to investigate the full Green's function in presence of 
the time-dependent electric field. 

The paper is organized as follows. In section we set up the general formalism (Dyson's 
equation), in section |IT| we specify to the plane wave basis and derive the exact expression 
for the Green's function which is evaluated numerically. In section |TV|, we give a formulation 
in the eigenstate basis of the static potential and discuss the electric field induced scattering 
rate, which again is derived from the formally exact result. A short conclusion is given in 
section |V]. 



II. GENERAL FORMULATION 
A. Dyson equation 

The starting point in each of our approaches is the Dyson equation for the Keldysh 
Green's function matrices. We start from a basis of eigenstates labeled with a and de- 
fine Green's functions |T3| iGlp{ti,t2) := {Tca{ti)c\{t2)) , iGlp{ti,t2) := {f Caih) cl{t2)) , 
iG>p{ti,t2) := (c„(ti)cj;(t2)), iG<p{ti,t2) := -{cl{t2)ca{ti)), where T (T) denote (anti)- 
chronological time ordering. Furthermore, the brackets () denote as usual an average over a 
quantum state of the system (not necessarily the stationary state). The Green's functions 
are written in a matrix block 

G^(ti,t2) -G<ih,t2)\ 

G>{UM) -G^{tiM))' ^ ^ 

for which the Dyson-integral equation has the same form as in equilibrium theory, 

G(ti, t2) = G°(ti, t2) + j I dtdt'G\h, t)^{t, t')G{t', ta), (2) 

where the matrix S is composed of S^, — E^, S^, — S*^, in analogy to Eq. (0). Here, G° 
denotes the unperturbed Greens' function, e.g. iG^p{ti,t2) := {Tca{ti)c^p{t2)) , where the 
electron creation (annihilation) operators c^^''(t) are given in the interaction picture which 
will be defined according to the splitting of the total time dependent Hamiltonian H{t). 
This splitting defines the way the perturbation theory is performed and will be explained 
below. 

Since we are interested in a Hamiltonian a part of which oscillates with frequency Vt, it is 
useful to perform a Fourier analysis according to time 'center of mass' and relative coordi- 
nates ("Wigner-coordinates"), namely T = {ti + t2)/2 and t = ti — t2. This decomposition 
is defined according to 

G{ti =T + t/2,t2=T - t/2) = — E / duje-'^^'e'^^^Giuj, N) (3) 



GitiM) = 
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and correspondingly for and S. The inverse transformation is 



G{u, N) = r dte'^' r ^^iii^e-*^^^G(ti = T + t/2, = T - t/2). (4) 

J-oo Jo 27C 

Of special interest is the component = which determines the average over the 'center-of- 
mass' coordinate T. In particular, in the case of an equilibrium situation (no electric field), 
all components G{u!, N) with 7^ vanish because the Green's function depends on the 
relative coordinate t = ti —t2 only. 

Inserting Eq. (0) into the Dyson equation Eq. (||), a straightforward calculation yields 

G{uj, N) = G%uj, AT) + ^ G° L + Ni) S L + N - - N2 

N1N2 ^ ^ ^ ^ ^ 

/ N - No \ 

xG(uj + ^-^n,N2). (5) 



d{QT) . 



In this paper, we concentrate on the non- interacting case where the perturbation is a one- 
particle operator and the self energy is (Keldysh) block-diagonal 

An impurity average effectively introduces interactions among the electrons and the self 
energy becomes different from Eq. @. However, as long as no impurity average is performed, 
one merely has to deal with a one- particle problem and the integral equation Eq. together 
with Eq. (I) exactly determines the Green's function G{ti,t2). 

Because of the linear relation G^ + G^ = G^ + G^, a rotation to tridiagonal form can 
be performed in Eq. (^). One of the resulting equations will be used later, namely the one 
for the retarded Green's function iG^p{ti,t2) := 6{ti — t2)(ca(ti)c^(t2) + c^(^2)ca(ti)), which 
reads, using S^(t,t') = V{t)6{t - t'), 

t2) = G°'^(ti, t2) + J dtG^^'^'Xh, t)\/(t)G^(t, t2). (7) 



B. Time-dependent field 

In all what follows we assume that the system is subject to a spatially homogeneous 
electric field which oscillates in time with a frequency Q and is polarized in direction e, 

E{t) = eEocos{nt). (8) 

The associated vector potential A^(t) = —{eEo/Q)sm{Qt) couples to the momentum of 
the i-th electron via Pi ^ Pi — A(xi) — A'^(t), where A(xi) is included here to indicate that 
there can be other contributions to the vector potential, e.g from a static magnetic field (we 
use units h = e = c = 1 throughout). The additional energy through the electric field is 

F 1 
H,{t) = -^sin(fit)5:ev, + A^eA$(t), A$(t) := ——E',sm\nt), (9) 
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where m* is the bandmass, m*Vi = Pi — A(xi), and Ne the number of electrons in the 
system. We exphcitely note that although the electric field E(t) is homogeneous, in an 
alternative gauge a corresponding scalar potential would be linear in the space coordinate, 
namely ~ excos(f2t), an 'oscillating slanted surface'. 

By gauge invariance, the current density operator is j(x, t) = jo(x) — ^A'^(t)p(x), where 
jo is the (paramagnetic) current density operator in the absence of the electric field, and 
p(x) the electronic density. 

In second quantization and in an basis of eigenstates {\a)} with eigenenergies the 
electric field gives rise to an additional Hamiltonian 

H^it) = ^ sm{nt) J2{»\v\P)c% + NAm. (10) 

a/3 

Here, v = (l/m*)e(p — A(x)) is the component of the electron velocity operator in direction 
e. Note that the coupling is effectively to the velocity and not to the current density operator 
since we assumed spatial homogeneity of the electric field. 

Care has to be taken when calculating Green's functions in a specific gauge, as has 



been pointed out by Bertoncini and Jauho [15]. Since only gauge-invariant quantities are 



meaningful for comparison with experiments, one has to formulate a gauge-invariant theory. 
In fact, a version G of the Green's function which is invariant under a gauge transformation 
of the electric field can be recovered from the Green's function G in any specific gauge. In 



real space, this transformation is as follows ITSl: if R = (x + x')/2 and r = x — x' denote 



center-of mass and relative coordinates in real space and t = ti —t2 and T = {ti + 12)/2 the 
Wigner time variables, 

G(r, t, R, T) = exp (iw(r, t, R, T)) G(r, t, R, T), (11) 

where 

w{r, t, R, T) := f^^^ dX [0(R + Ar, T + At) - rA"(R + Ar, T + At)] , (12) 

J-l/2 



and and A"^ denote the static and the vector potential, respectively. Eq. ([11]) can be 
proofed directly by performing a gauge transformation A^ A*^ -|- Vx, — > — dtX and 
using the fact that the field operators in real space transform as \l/ ^ exp[ix)'^ (note that 
we have set e = ^ = c = 1). The factor exp{iw) becomes particularity simple in our case of 
a spacially homogeneous electric field defined through = and A'^(t) = —{eEo/Q) sin (fit). 
If we start in an arbitrary basis {|a)} with wave functions \I'q,(x), the gauge invariant object 
is 

G(r, t, R, T) = 5: v[/^(R + r/2)v[/;(R - r/2) exp (^^(r, t, R, T)) G*^(t, T), (13) 

a/3 

where G'*^(t,T) is calculated in a specific gauge and basis \a). 

In the rest of this paper, we will always refer to a specific gauge and refer to Eq. (p!3D 
for the final, gauge-invariant form. The importance of the phase factor exp{iw) depends 
on the physical situation and has been discussed, e.g., in p!5|,p!6|. For convenience, we split 



off another phase factor in Eq. (|T^). The last term in the Hamiltonian Eq. ([T0|) which is 
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quadratic in Eq couples to the total particle number operator N and thus gives rise to a 
phase factor in the Green's function. The latter can be written as 



G'*(ti,t2) = e 



N 



t2 



dt'A(^{t'). 



(14) 



where is the particle number and G the Green's function for the same Hamiltonian 
Eq. (|l^) but without the last term which, as can been directly seen from Eq. (10), shifts all 
energy levels by the same, but time dependent amount A$(t). 



III. PERTURBATION THEORY IN THE IMPURITY POTENTIAL 

In the preceding section we have derived the general equation Eq. (^ which determines 
the Green's function of the system. As mentioned above, the splitting of the total Hamilto- 
nian into perturbation part and unperturbed part defines the unperturbed Green's function 

which is at the base of any perturbation theory. In the following, we define the model 
Hamiltonian. In this section, we start from an eigenstate basis of plane waves, i.e. |a) = |k) 
where (x|k) = = (l/L*^/^) exp(— ikx) and L"^ is the system volume (L ^ cxd in the 

thermodynamic limit). We furthermore assume that there is no magnetic field, a non-zero 
magnetic field can be included in the basis of exact scattering states treated in the next 
section. In the plane wave basis, the velocity matrix element is diagonal, 

{a\v\P) = Sa,p{a\v\a) = 5apVa, (15) 

namely Va = = ek/m*. Although Eq. (|I5|) is exact only for plane waves, we keep the 
notation as general as possible and use greek letters labeling the eigenstates. The following 
results are then also valid for a general eigenstate basis with the approximation that non- 
diagonal elements of the velocity {a\v\(3) are zero. Because of Eq. ([ISD , the Hamiltonian 
(system -|- electric field) is 



Ho{t) = J2 ^c,{t)cic^; Bait) =ea + -^Va sin(nt), (16) 



where Ea is the energy of state a. 

The static potential V^(x) now gives rise to an additional part to the total Hamiltonian 
which reads 

Hit) = Hoit) + V, V = Y.V^pcicp. (17) 

a/3 

At this stage, V is not yet specified further. Depending on the physical situation, it describes 
a single impurity, a distribution of random scatterers, a double barrier etc. Note that 



according to the above, in Eq. (17) a denotes a plane wave while in the corresponding 
'diagonal approximation' case for arbitrary basis a the potential term V does not appear 
and simply H{t) = HQ(t). 
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A. Dyson equation in arbitrary strong electric fields 

We define an interaction picture with respect to the unperturbed part Hq of the Hamil- 
tonian according to 

c„(t) = U^o{t,to)cMt,to), Uo{t,to) = Texp (^-i j' dt' Ho{t')^ . (18) 

The unperturbed time-ordered Green's function then is 

iGl;;{hM) = {Tc^{h)cl{h)) = 

= e-'^^^'"''"^'^ [Bit, - t,){l - - d{h - 6^^, U = (cicj (19) 

and correspondingly for the other Green's functions. The average () gives rise to a distribu- 
tion function which not necessarily coincides with the Fermi distribution. Note that the 
value of the fixed time to in Eq. ([T8|) is not relevant for G. The Green's function Eq. ( p!9[) 
can be related to the free {V = 0) Green's function without electric field {Eq = 0), which in 
the following is called G^"''^'^. Namely, one has for the time-ordered part 

G^f,{hM) = e^9.lcosm)-cosint,)]GlJree^t, - 1,)5^,, .= ^ (20) 

and correspondingly for the other Keldysh-block components. The effect of the electric field 
thus is incorporated in the phase factor in front of G'^'^^^'^ in Eq. (pO]) . Note that this phase 
factor, in contrast to the 'global' phase $(^1,^2), Eq. (0), depends on the quantum number 
a which will be important when scattering, i.e. transitions between different states a —>■ (3, 
is considered. 

To get into contact with the Fourier representation above, we use a decomposition into 
Bessel functions according to exp[iz cos{rit)] = X^n ^"'^n(^) 6xp(mnt) (the sum is over all 
integers n). We obtain 

Glp{u,N) = z''Y.U9a)JN^n{-9a)Gi'''''{u;+[n-N/2]Q)6o^p (21) 

n 

for the Keldysh Green's function matrix. Note that G*/^*^*^ depends on one Fourier variable 
only since the corresponding G-^^''^'^{ti — ^2) is a function of the time difference only. We 
note again that the eigenstate basis in this section consists of plane waves |a) = |k) with 
momentum k. Only in the 'diagonal approximation' model, where one neglects the potential 
term V in Eq. (p!7|), the results of this section can be transfered to an arbitrary eigenstate 
basis. Furthermore, in this section G° denotes the Green's function belonging to the time- 
dependent Hamiltonian Ho{t), Eq. ([T6|). It contains the electric field exactly to all orders of 
the coupling parameter ga, Eq. (|20|) , which appears in the argument of the Bessel functions. 

We now proceed to set up a systematic approximation scheme that allows to calculate 
the Fourier components G^fsiuj, N) of the Green's function belonging to the time- dependent 
Hamiltonian H{t), Eq. (0). In particular, the effect of the time dependent field will be fully 
incorporated through the unperturbed Green's functions G^, Eq. (|2T|) . 

Since the perturbation potential V is time-independent, the selfenergy in Eq. (|^) is 
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^u,N) = S-6n,o, 5= ^j. (22) 

Thus, 

G{u, N) = G\iu, N) + J2G^(^^ + iVi) SG (^u + - A^i) . (23) 

One immediately sees the fundamental difficulty in Eq. (PB|): Even for obtaining only the 
Fourier component G{u},N = 0), corresponding to an average of the 'center-of-mass' time 
coordinate, the A^i-sum couples the different A^-components in Eq. (^). 

An obvious approximation would be to neglect all A^ 7^ components at all, i.e. to 
neglect the center-of-mass time coordinate. Indeed, approximations of this kind have been 
used recently (we call it 'fast' approximation in the following), namely in the so-called 'non- 
adiabatic regime' of resonant tunneling where the frequency Q of the time dependent field 
is much larger than the inverse tunneling time. Another example is the neglect of higher 
order Fourier components in the calculation of Green's functions for Wannier Stark systems 
rT|. There, the approximations were argued to work well for frequencies Q larger than the 



bandwidth of the systems (the letter is proportional to the hopping matrix element of the 
corresponding tight-binding Hamiltonian) . The extreme advantage of those approximations 
is that in a relatively simple manner the effect of the time dependent field can be incorporated 
by essentially calculating expressions where the free Green's function G^^'^^iu) is replaced 
byGO(u;,Ar = 0), Eq. (|21|). 

In the general situation of an (arbitrary) disorder potential or an arbitrary eigenstate 
basis, it is not clear if an analogy to the fast approximation exists. Furthermore, the validity 
of the latter has not been investigated systematically to our knowledge. In the following, we 
develop a scheme for such an investigation. The formalism set up below is completely general 
and can be applied to a large range of other situations, too. This opens the possibility to 
go beyond the 'fast approximation' in a general scattering situation. 



B. Fast approximation and its first corrections 

We start calculating the A^ = 0, ±1 components of G{uj,N), Eq. (|23|) , i.e. the 'fast' ap- 
proximation which neglects the 'center-of-mass' time coordinate T, and its ffist corrections. 
Introducing the convenient notation 

^S:=G(°)(^ + ^l],n), (24) 

we obtain from the terms A^^i = 0, ±1 in Eq. (^) the approximate solution 

9l,n ~ g'ln + 9^-nfiS9l,n + 9l+l-n,lS 9l+l,n-l + 9-l+l-n,-lS 9l-l,n+l- (25) 

Here, all quantities are (Keldysh block-) matrices. The problem of determining the (time- 
average) component go^ = G{uj, A^ = 0) is pure algebra and one obtains 

G{u;,N = 0)^{l-~g'S)-'~g' 

~9' = <o + 9lS{l - g%Sr'gl_, + g\_,S{l - g^^oSr'g^^. (26) 



8 



Equivalently, other components than = might be calculated. Note that Eq. p6| ) has 
the standard form of the formal solution of the time-independent Dyson equation for a self 
energy S (which for our case of a static potential V is given by V itself), 

G{uj) = (1 - G^'''''{uj)S)-^G^''^'{uj). (27) 

However, the 'unperturbed' Green's function g^Q is substituted by the more complicated 
expression In fact, having g^Q instead of in Eq. (^) is exactly the 'fast approximation', 
namely 

G{uj, N = 0)^(1- G%u, N = 0)S)-^G\uj, N = 0). (28) 

Eq. (^) is the zeroth order approximation for Eq. (^) with respect to the Fourier index 
A^. In particular, the dependence on the 'center-of-mass' time coordinate T, Eq. (|^) is 
completely neglected here. 



C. Exact solution 

The perturbative calculation of the = 0, ±l-components of G{uj, N) above suggests 
that there is a general solution of the Dyson equation type which holds to all orders. Eq. (|23D 
can be solved formally to all orders in the self energy and for all components N. The proce- 
dure for this is very similar to the derivation of the Floquet-matrix in the exact eigenstate 
basis below. Here, we derive the exact result in the plane wave basis. It is the generalization 
of the Dyson equation solution for a time-independent scattering problem to the case of an 
additional time-dependent electric field and will be the starting point for the subsequent 
numerical truncation scheme. 

We first rewrite Eq. ( PB| ) as 

9ln = din + H 9^+n',n'+nS9l+n+n'-n', (29) 
n' 

where we shifted the n'-summation index. We introduce a matrix 7 with elements 7^5 = 
gin = gr+s,s-r, together with a matrix 7° with the corresponding elements of g^. We perform 
an index transformation 

r = (/-n)/2, s = (/ + n)/2, (30) 

which yields 

7.. = 7°. + E7°n'^7n'.. (31) 

n' 

Therefore, in the space of the (r, s)-indices, the solution of the matrix equation becomes 

7 = (l-7°5)-Y. (32) 
Explicitely, the matrix 7 has the form 
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X 7°. (33) 



I - 

7= ••• -^-i-i-^ I-Qq^S -g'^^S ... 

-9L2S -gUs 1-9%S ... 

\... ....J 

In particular, the matrix (1 — 7°S') is fully occupied. 

Since 7 is just another representation of the original matrix g, ■jrs = gr+s,s-r, Eq- (|33D 
solves the Dyson equation Eq. ( p3D exactly. Note that the elements of the matrix 7 are 



(Keldysh) matrices themselves. This tensor structure is analogous to the formal solution 
in the exact eigenstate basis below, here we started from a plane wave eigenstate basis |k). 
While 7° is diagonal in this basis, the self energy S in general is not diagonal and depends 
on k and k' for a given static potential. 

In the usual time-independent case, solving the Dyson equation Eq. ( pTf ) for an arbitrary 
potential requires the inversion of a matrix of dimension D, where D is the dimension of 
the eigenstate basis (which is finite in any numerical calculation). In the present time- 
dependent case, the effective dimension of the matrix 7 is D times Uph, where Uph is the 
number of 'photoblocks' in (1 — 7''5'), e.g. Uph = 3 for the three blocks per row in Eq. (p3D. 
Since there are infinitely many Fourier components of the Green's function, in the exact 
solution Eq. (|33|), Uph = 00, and 7 is of infinite dimension even if the dimension D of the 
original Hilbert space is finite. It turns out, however, that it is sufficient to consider only a 
relatively small number of these blocks in practical calculations. In principle, the numerical 
effort increases rapidly with the number riph. Fortunately, the examples below show that 
convergence is reached quickly by truncating 7 at relatively small Uph. We also stress that 
this 'truncation method' is to infinite order in the static potential V and the electric field 
Eq. It is perturbative in the higher Fourier components (large n) for the center-of-mass 
time coordinate T in gf^ in the off-diagonal in Eq. (p3[). 

Since Eq. (|3^) is the exact solution of Eq. (|23|) , one easily recovers the above approx- 
imations: The 'fast approximation' Eq. (^) is obtained by cutting the matrix (1 — 7°S') 
such that only the (r = 0, s = 0), i.e. the 'central' element is retained. In (r, s) space the 
(infinite) matrix becomes a number and the inversion is trivial. The next approximation 
Eq. (p6D that retains only the = 0, ±1 components can be verified similarly by inversion 
of the corresponding 3x3 matrix in (r, s)-space. 

Eq. (|33D offers the possibility of a systematic investigation of the 'truncation method', 
in particular of approximations like Eq. (|26|) or Eq. (|28|) . 

D. Numerical examples 

In the following, we apply the formalism to calculate the time- averaged Green's function 
G{u!, N = 0) numerically, using Eq. (|33D in two different potentials V. 



1. Constant potential 

In the first model, we use a constant self energy V = const. Physically, this is a trivial, 
nevertheless very instructive situation since it offers comparison to an exact solution. That 
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is, the constant shift of the energy scale by V corresponds to the exact result 



G,iu,N) = Glico-V,N) (34) 

which has to be reproduced by inverting the full Dyson equation, Eq. (P5|). Here, G° is the 
unperturbed (V = 0) Green's function defined according to Eq. ( pi]) . 

The result of the numerical inversion is shown for the N = 0-component of the corre- 
sponding spectral function component which we defined from the Green's function as 

Ak{uJ,N = 0) := --lmG§{uj,N = 0), (35) 
vr 

through the retarded Green's function. The numerical calculation here is regarded as a 
check for the consistency and convergence of the formalism as such. Therefore, for reasons 
of simplicity we did not calculate the full 'physical' spectral function Ak{uj) JlBl which 



involves the additional phase factor w, Eq. (^3|). Note that because of Eq. (^, we only have 
to consider retarded functions in the Dyson equation. In particular, in our plane wave basis 
the elements 1 — G^^S = 1 — G^^ are (complex) numbers and not matrices since we are 
working in Fourier space and not in real space here. 

The results of the numerical inversion of Eq. ( |33D are shown in Fig. (|1|) and Fig. (^. 
We fixed the wavevector k = kp, energies are measured in units of the Fermi energy Ep = 
{kkpy /2m* . The dimensionless parameter a regulating the coupling to the electric field of 
strength Eq is defined as 

a = EokF/{m*n^). (36) 

For a potential V = 0.22Ef and strong coupling, a = 2, the central peak of the spectral 
function is strongly suppressed: its main weight comes from the Bessel functions Jq, on the 
other hand, the product ak is close to the first zero of the latter. Side-peaks appear at values 
of the energy u shifted by multiples nQ. The exact result Eq. (p4D is excellently reproduced 
by inversion of the 11 x 11 'photo'block matrix Eq. (|33|), while the 'fast approximation' 
Eq. ( pHD gives a totally wrong result. In contrast to this, already the inversion of the 5x5 
matrix qualitatively reproduces the exact result. The situation is less dramatic for smaller 
coupling a = 0.5 and smaller potential V = O-OSEp. There, the 'fast approximation' Eq. ( |28D 
is qualitatively roughly correct whereas the exact result is already obtained from the 5x5 
matrix inverse. 



2. Resonant tunneling 

In the second example, we investigated a physically non-trivial situation, namely the 
transmission of an electron through a double-barrier structure. The model is a one- 
dimensional lattice with j^ax sites and a potential 

S = V = \x,)V,{Xi\ + \X2)V2{X2\. (37) 

For Vi,V2 > 0, this potential has a resonant structure in its transmission coefficient and 
can be regarded as a simple version of a model simulating a quantum well double barrier. 
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The Dyson equation Eq. (^) has to be transformed to real space to obtain the matrix 
elements {x\G{uj, N) \x') of the full Green's function in real space. Changing from the original 
plane-wave basis to real space in particular involves a Fourier transformation in the Green's 
function 

G\x - uj,N) = ^ f dke'^^''-'''^Gl{u, N). (38) 

ZTT J 

In fact, the integral Eq. (38) is tedious and we found an analytical form for the retarded 



Green's function only for small couplings xq = Eq/ 171*0,"^, namely 

G'^'-{x-x';u,N) = [sgnix-x')^fJ2Jni^o^/^)■JN-n{-Xo^/^)^e'^'\^-^\ 

z:=uj+{n-N/2)n, 2xo<\x-x'\ (39) 

For smaller |x — we resorted to a numerical evaluation of the integral. We defined the 
coefficient 



T{u;) 



G''{x,x',uj,N = 0) 



G^'^ix- x',LU,N = 0) 



(40) 



from the time-averaged = 0-components of the Green's functions which is independent 
of X, x' for X < Xi and x' > X2, i.e. transmission of an incoming electron (from left) through 
both potentials at sites xi and X2. Strictly speaking, since T{uj) is defined by the product 
of the time-averaged G's and not the time-average of the product GG*, Tioj) deviates from 
the (physical) time-averaged transmission coefficient. Here as above we restrict ourselves, 
however, to the simpler quantity Eq. (^) to clarify the convergence of the numerical method 
to determine G. We chose jmax = 12 and a distance X2 — Xi = 2a, where a is the lattice 
constant serving as length scale and Ea = fi /2m* a? for the energy scale. The number nph 
of 'photo-blocks' in the matrix (1 — 7''5') again determines the total size of the matrix to 
be inverted. Each block itself is a jmax x Jmax^niatrix, the total dimension of the matrix 
therefore being jmax x riph. 

Fig. d^) shows the coefficient T{uj) for three different sizes of the matrix Eq. (pSj). One 
reckognizes that already for = 5 convergence is reached practically for all values of 
the energy u. The side-peaks to the left and to the right of the transmission maximum 
correspond to the first 'side-bands' ±hQ for tunneling through the double barrier and have 
been found previously by Inarrea et al. [0 . 

We conclude that the numerical inversion of Eq. (|33D shows the extremely good conver- 
gence of the results from truncating Eq. (^), in contrast to the approximation Eq. (^B|). 
This in particular indicates that one has to be very careful when using approximations where 
Fourier components of the Green's function are neglected. We rather believe that our cal- 
culation shows that in general one has to consider at least a certain number of components 
G{uj, N), N 0, from the very beginning. In the original representation of G{ti, ^2), Eq. (^, 
this means that not only the relative time coordinate ti — ^2 but also the 'center-of-mass' 
time coordinate T is important. 
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IV. TIME DEPENDENT FIELD AS PERTURBATION 



The basic advantage of the above approach is its exactness to all orders in the electric field 
which a priori is incorporated in the free Green's function G^. The exact result Eq. ( |33[ ) 
implies that one solves the static scattering problem and the effect of the electric field 
simultaneously. By inverting the truncated matrix one sums up the static potential V to 
infinite order. Still, higher Fourier component terms Sg'^g proportional to V in the off- 



diagonal of Eq. (|33D are neglected. On the other hand, there are situations where one wishes 
to make no approximations concerning the static potential, in particular in systems where 
localization plays a role. Furthermore, in cases where the static potential is complicated, 
practical calculations become tedious in the plane wave basis. In fact, if the static scattering 
problem is already solved (say by numerical diagonalization), the basis of the static potential 
eigenstates provides the ideal starting point. 

In the following, we therefore take a different point of view and consider the electric 
field as a perturbation while starting from the exact basis of eigenstates of the system, 
i.e. the scattering states of the static potential V. In presence of an additional, static 
magnetic field B, these are the scattering states of the static potential in presence of B. 
As above, a formally exact expression for the Green's function is derived which is closely 
related to the Floquet state formalism . As an application, we derive an exact expression 



for the dephasing time due to the time-dependent electric field expressed in terms of the 
AG-conductivity of the system in absence of the electric field. This expression in particular 
is valid near or at a metal-insulator transition. 
We start from a Hamiltonian 

Hit) = Ho + M{t) = Y: e^cic, + ^ M,^(t)4c^ 

a a/3 
771 

M„^(t) = Maf3 sm{Qt) = -^{a\v\(3) sin(fit). (41) 

The unperturbed part Hq is now time-independent and describes the non-interacting system 
in the exact eigenstate basis, labeled with the index a. An additional advantage is that the 
unperturbed quantities are now defined in an equilibrium system in contrast to the case in 
the preceding section. There, Ho{t) was time dependent and thus even the unperturbed 
system not in equilibrium. 

The matrix element Mq^ describes transitions from states a to states (3 and couples to 
the velocity matrix element in polarization direction of the electric field. In Fourier 

representation, the free Green's function G^ belonging to Hq is diagonal and depends on 
the time difference only, while the self energy T,{t,t') = 6{t — t') sin(f2t)E has components 
N = ±1, 

G%iu,N) = 5^,oSapG'iiuj) 

E^fsiu, N) = S„^[5;v,i - 5^,-1], ^afs = ^W^a/3 ^ U ^ ] . (42) 



Using the general Eq. (13) one obtains 



Giuj,N) = G\iu)6N,o + G'' 



uj L 

2 ) 



af3 



G{lu + ^,N-1]-g(lu-^,N + 1 



2 / V 2 



(43) 
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For convenience, we again use the notation Eq. (^) gi^n = G{uj + lQ/2,n), furthermore 
gf = G'^(uj + m/2), so that Eq. (||) becomes 

9l,n = ^nfiQ^l + - gi-l.n+l]- (44) 

We note that again all quantities are Keldysh block matrices, each block being a matrix in the 
indices that label the exaci eigenstates of Hq (static scattering potential plus magnetic 
field). The free Green's function G° in this chapter thus refers to the time-independent part 
Hq only and is diagonal in the eigenstate basis. 



A. Exact formal solution 



We use the index transformation Eq. (|30|) r = {I — n)/2, s = {l + n)/2,r, s = 0, ±1, ±2, ... 
which allows for a formal solution of Eq. (^i]) and shows the relation to the Floquet state 
theory. Again, we define the matrices 



2r 



9r+s,s 



T-0 



(45) 



Note that in contrast to the previous section, F and g now refer to the exact eigenstate 
basis. The requirement that both r and s are integer numbers in fact restricts the value of 
/: Even n = s — r requires even / = s + r and odd n requires odd /. From the definition of 
9in = G{uj + l/2Q,n) one reckognizes, however, that the index / is redundant in the sense 
that one value of I is sufficient to obtain the full information of the function G{uj,n). This 



kind of redundance is well-known from the Floquet state formalism [|r^ to which the analogy 
becomes evident in the following. 

From Eq. (0), we obtain a corresponding equation in the r, s— space, using gi±i^nTi = 

F. 



F^g + F2j,S [Fr+l,s ~ Tr-l.s 



After introducing 



this can be written in the matrix form 



(5r-l,r')^! 



F = F° + FVF = [(F 



a 



(46) 



(47) 



(48) 



Eq. ( ^81 ) is the formal solution of the problem to determine the Green's function for the 
Hamiltonian Eq. (^Tj). The structure of the matrix F is as follows: Since F*^ is diagonal in r 
and s, one has (F°)~^ = 5rs{g^=2r,n=Q)^^ = drsG^iu + rVt)^^. Therefore, the matrix F is the 
inverse of a tridiagonal matrix. 



s G\Loy^ -s 







-E 

G\uj + 2n)-^ -S 
s ... / 



(49) 
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i.e. an infinite tridiagonal block matrix in the space of 'pfioton numbers'. Note tliat tlie 
elements of 7 are (Keldysli) matrices themselves and in general infinite-dimensional as is the 
Hilbert space belonging to the Hamiltonian Hq. 

To give an illustrative example, we consider a two-level system with energies Sa and e^, 
coupled to a time-dependent field via the matrix element Map = Mg^ = 26, M^a = Mpp = 
0. Taking the limit ^ in Eq. (^Of), we obtain (notice that the quantity F is just a 
representation of the Green's function G itself): 



^^(a; ^ 0) 



\ 



-1 



Q + Ea 

n + Ep 

ib 

ib 



-ib 



-ib 





-ib 



V 



(50) 



In the case of a coupling proportional to cos{Qt) instead of sin(f2t), one easily verifies that 
±ib is replaced by b in Eq. (^O]). One then has exact coincidence with the formal solution 
of the time-dependent two-level problem first given by Shirley |]l7 in his seminal paper on 
the Floquet state theory. Here, we have reproduced those results in the language of Green's 
functions starting from the Dyson equation. The Floquet Hamiltonian Ti is obtained as the 
limit Ti. = —{G^)~^{uj 0) in the representation Eq. (^5]) without directly refering to an 
extended Hilbert space. 



B. Dyson equation in second order in the electric field, electron scattering rate 



We consider Eq. (|4^) to derive the first and second order expression for the Green's 
function G. The first order response to the perturbation M is obtained from qq^i, which 
is an example of how to reproduce results from standard equilibrium theory (Kubo's linear 
response expression) from Keldysh Green's functions. To next (quadratic) order in E, we 
obtain the = (time averaged component) G{uj, N = 0) from 



fi'0,0 

which can be written as 

G{uj,N = 0) = 
S{uj) 



(51) 



(G 



0\-l, 



S{UJ] 



-1 



^WG\uj + Q)W + {Q 



(52) 



Specializing to the retarded functions, Eq. (|^), we find G^{uj,N = 0) = [{G^'^) 
S^{uj)]-\ where S^{uj) := (l/4)MG°'^(cj + fi)M + {n ^ -^]). 

The determination of the Green's function Eq. ( |52D requires a full matrix inversion since 
S in general is a non-diagonal matrix. The diagonal part of the imaginary part of the 
retarded self energy ImS^a^ however, has adirect physical meaning of an inverse lifetime 
of the state a. It describes the scattering rate for scattering out of state a to any other 
eigenstate of the system due to the perturbation. This quantity can be expressed exactly by 
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the real part lie (Ti{e = ep,^) of the AC conductivity as a function of the Fermi energy ep 
and frequency Q in absence of the time-dependent perturbation Eo(t). In the hmit of zero 
temperature T = 0, the derivation given in App. |^ yields 

r^-r-^.r:', 4^^^^^^.^^ + O f^V . T ^ 0. (53) 



2 (hnypie) 

where we reinstalled the h and p{e) denotes the density of states. This equation relates a 
physical quantity (the scattering time t±) of the non-equilibrium system (static potential 
+ electric field) to a linear transport-quantity of the equilibrium system (no electric field). 
The validity of this relation is independent of the eigenstate basis and in particular non- 
perturbative in the impurity potential. Basically, the reason why both quantities can be 
related to each other is that both essentially are given by the square of the matrix element 
of the velocity operator. It is also obvious why the (Fermi) energy e is shifted by ±hQ/2 
in the argument of the conductivity. The latter is due to particle hole excitations that lie 
symmetrically around the Fermi level e, whereas t± describes electron scattering from e to 
e±h^l. Furthermore, the expression for T^]y, Eq. (0), is valid for a single electron ('empty 
band'), i.e. no Pauli block by the other electrons. If the Fermi sea of the other electrons is 
considered, only scattering through photon absorption (+) is possible at zero temperature, 
and TJ^J^r = t_^_^. 

The physical meaning of Eq. (^) becomes particularily transparent by recalling that the 
power Pq absorbed by electrons subject to a homogeneous electric field Eq is given by 

Pn{e) = {l/2)E^o^eai{e)-L'', (54) 

where L*^ is the system volume. For frequencies fl such that hQ <^ e, one can set 6±hQ ^ e, 
and Eq. (^) can be written 



hn A{e) + 



s)=p{e){L''hn)T-\e), (55) 



where A(e) := {L'^p{e))~^ is the mean level spacing at the energy e. The number of photons 
absorbed per time, Pn{£)/hQ, is given by the scattering rate, multiplied with the density of 
electronic states p{e) on the scale of the photon density of states {L'^hQ)~^. 



C. Microwave radiation in the integer quantum Hall efFect (IQHE) 

As mentioned in the introduction, in the integer quantum Hall regime, microwave induced 
deviations AR^x, AR^y of the resistance tensor have been observed experimentally ^j. One 
main finding was that these deviations are proportional to the square root of the incident 
microwave power, i.e. linear in the electric field amplitude Eq. 

We calculated the microwave induced scattering rate Eq. (|53D and found that this linear 
dependence cannot be explained in a simple scaling picture through the interplay of local- 
ization length and inelastic length. In fact, in the past the standard scaling picture of the 
IQHE [l^ has been widely used to describe the scaling of the temperature (T)-dependent 
slope (dpxy/dB)"^"'^ and the inverse of the half- width of the pxx{B)-peaks |T^-^, i.e. 
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temperature-dependent deviations AR^x/xy The main idea is that the IQHE as a metal- 
insulator transition can be described in terms of a scaling function, depending on the ratio 
^/Lini where ^ is the localization length and 

= {Drf/' (56) 

is a T-dependent inelastic length scale due to electron-electron and/or electron-phonon 
scattering at a rate r^^. D is a diffusion constant. 

Although it is tempting to use this picture under inclusion of the additional inelastic 
scattering process due to the microwave, we found that such an explanation definitely is not 
appropriate. Our argument is as follows. At low temperatures, electron-electron scattering 
is believed to be dominant; in presence of microwave radiation, however, the scattering 
described by the rate Eq. (^31) also should contribute. If the scattering due to the microwave 
is the main inelastic process, t]^Iw ^ ^ where stands for the rate due to electron- 
electron, and r^pj/ for the scattering by the time- dependent field, Eq. (|53|) . One then has 
Lin ~ {DtmwY^'^ Eq. (|56D, and the scaling variable ~ Eq becomes proportional to 

the amplitude Eq of the microwave field which would hint to the linear Eq- dependence of 
ARxx/xy On the other hand, the electron-electron scattering is dominant for ^ TkIw 
Estimating the effect by the microwave by adding the scattering rates, = + Tj^^, 
one finds in this case 



r-' = r;^\l + 0{E'q)) . (57) 

In particular this would mean that the (small) correction to the scaling variable due to the 
microwave radiation is quadratic in the amplitude Eq, and so should be any deviation to the 
resistance Rxx/xy- 

We used the parameters of the experiment 0, Eq. (|53|) , and a previous calculation of 
the rate t~^^ p2[ to compare both quantities. It turns out, that for temperatures T ^ 0.5K, 



fl = iOGhz, Eq ?a lO'^V/m, the rate t^vf orders of magnitudes smaller than the electron- 
electron scattering rate t~J^ = 'ykBT/h |]22|, where 7 is a factor of order unity. Numerically, 
we find values approximately as r^"^ ^ lO^^s^^ and T^j\y ~ 10^s~^ for an estimated (small) 
conductivity axx ~ I0~'^e'^/h. Even for a conductivity axx close to the 'ideal' value 0.5e'^/h, 
the rate t^\v stays an order of magnitude behind the electron- electron scattering rate, and 
Eq. (|57|) should apply. This, in turn, is in contrast to the experimental finding AR ~ Eq. 
From this negative result we conclude that an interpretation of the experimental data 
is not possible in the above picture, simply using scaling arguments in combination with 
inelastic scattering rates. We believe that one has to go beyond an explanation that only 
makes use of scattering rates, i.e. the diagonal part of the self energy. Rather, one should 
start from the full Green's function Eq. (^) and investigate directly how quantities like the 
localization length are affected by the time dependent perturbation. 



V. CONCLUSION 

In this paper, we developed a method to calculate the Green's function of noninteracting 
electrons in static potentials in presence of time-dependent electric fields. The formalism 
was given in terms of Keldysh matrices (which are block-diagonal for static scattering). 
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Although not discussed here, in principal it allows inclusion of interaction processes in form 
of a suitable self energy for electron-electron or clcctron-phonon interactions. 

In the first approach, by a double Fourier expansion a systematic perturbation was 
developed which is exact in the electric field and sums up the static potential to infinite 
order. Numerical examples for a constant and a simple double barrier potential showed 
that a few Fourier components already are sTifRcicnt to obtain convergence. Wc furthermore 
showed that this 'truncation method' is a systematic expansion around a non-adiabatic limit 
of high frequencies Q of the electric field. 

A second approach was given in the basis of the exact scattering states. There, the 
truncation method is a perturbation theory in the electric field which was explicitly shown 
for the self energy to second order. The latter was used to derive a relation between the 
rate for (inelastic) scattering induced by the time dependent field, and the AC conductivity 
in absence of the field. Application to microwave experiments in the integer Quantum Hall 
effect suggested, however, that rather than the microwave scattering rate the full Green's 
function in the basis of the scattering states has to be investigated. 

Although in both methods exact solutions were derived, the numerical applicability sug- 
gests that the first approach is appropriate for relatively simple static potentials and strong 
electric fields. The second should be used for arbitrary (e.g. random) potentials but rather 
weak electric fields. Although the second approach has not yet been tested numerically, we 
propose to apply it to the localization-delocalization problem. 

Discussions with B. Kramer, A. Kawabata, G. Platero, and L. Schweitzer are acknowl- 
edged. This work has been supported by the EU STF 9 program and the CREST program 
of the Japan Science and Technology Corporation. 
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APPENDIX A: DERIVATION OF THE SCATTERING RATE r^|^ 

Here, we derive the expression Eq. (^) for the inverse lifetime of a state a due to 
the time dependent electric field. First, the unperturbed Green's function is G^p{uj) = 
Sa/si^^ — £a + iO'^)~^. Note that energies are counted from zero and not from a chemical 
potential here. The coupling matrix element to the electric field is Majj = {Eo/Q){a\v\f3) , v 
is the component of the velocity operator in direction of the electric field polarization. This 
yields 




Y^{a\vh){^\v\P)5{cu-e, + n) + {n 



7 



(Al) 
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where we defined the components corresponding to photon absorption {+Q ) and emission 
(-^]) in Eq. 

We now recall the expression for the real part of the longitudinal AC conductivity in 
linear response to an arbitrary homogeneous electric field of frequency Q, namely in the 
zero-temperature limit 

Re aiiep, f^) = ^ E I {a\vm'S{eF + n/2 - ep)S{eF - n/2 - e^) + O (—)' . (A2) 

^ af3 ^^^^ 

We explicitely indicate the dependence of ai on the Fermi energy, furthermore v is the 
component of the velocity operator in one fixed direction, say Vx for cr^x- We assume equality 
of all diagonal components of the conductivity tensor and the limit of zero temperature; note 
that the expression explicitely implies symmetry in Q as it must be since Re crj(ep', f2) = 
Reai^Ep, — fi). 

The rate for scattering out of a state a with energy e = Ea due to photon absorption (+) 
or emission (— ) is given by the diagonal part of the self energy, 

rMW = r^' + r:\ ri\e) = -2lmS^^ie = e^). (A3) 

The relation between and the self energy is as in the equilibrium case and can be 
directly verified by calculating the transition probability from a state a to a state f3 due to 
the time-dependent potential M{t) as in the derivation of Fermi's 'golden rule'. 

We assume that the diagonal elements of depend on the quantum number a only 
through the corresponding eigenenergy Ea- In particular, for a (spin-polarized) disorder 
broadened Landau band discussed below, there is a one-to-one correspondence of eigenstates 
and eigenenergies. From Eq. (|A2|) , one finds 

Reai{eF±n/2,n) = 7rp(e^) J] | ^^(e^ ± - e^) , (A4) 

^ £a=£F 

where p(e) = {l/L^)J2a^{^ ^ ^a) is the density of states. Combining Eq. (|Al|), Eq. ([A3[), 
and Eq. (|A^), one obtains Eq. (E 
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FIGURES 



FIG. 1. The = spectral function component Eq. ( |35D in units of the inverse Fermi energy 
Ep^ at the Fermi wave vector k = kp- The model system describes a constant potential V and a 
time-dependent electric field of frequency fi. The dimensionless parameter a regulates the coupling 
to the electric field, energies are in units of Ep. The exact solution is shown together with the 
result from the 'fast approximation', Eq. (^), and the result from the numerical inversion of 



Eq. (33). Already for Uph = 11, practically for all energies coincidence is reached with the exact 



result Eq. (34) 



FIG. 2. Same as in Fig. (|T|) but for weaker coupling a and smaller potential V. Coincidence is 
reached for an even smaller matrix Eq. (^) in this case. 



FIG. 3. Transmission coefficient Eq. (|4(3| ) through a one-dimensional double barrier model in 
a time-dependent electric field with frequency O and coupling parameter a. Again, the Green's 
function has been calculated numerically for different numbers Uph of 'photoblocks' in order to check 
the convergence of the results from inverting the truncated matrices Eq. ( ^3| ) . The system is defined 
by a Id tight-binding model with 12 sites. The distance of the two potentials of strength V = 20Ea 
defining the barrier is 2a with a the lattice constant; energies are in units of Ea ^ f)} j^rrC o? . Site 
peaks in the transmission coefficient occur a energies shifted by ±^$7 from the central transmission 
peak. 
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